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ATMOSPHERIC EFFECTS ON CO^ LASER PROPAGATION 

I. INTRODUCTION 

The Marshall Space Flight Center (MSFC) has been involved in the 
development of CO 2 Laser Doppler Velocimeter Systems for atmospheric mea- 
surement since the late 1960's. The shorter range continuous wave (CW) 
system development has led to a variety of applications, including measurement 
of the following: aircraft wake vortices [1] , dust devils [2] , smokestack 
emissions [3] , and wind profiles [4] . The development of a pulsed system has 
resulted in measurement programs for clear air turbulence (CAT) [5J and 
severe storm research [6]. 

Accompanying the development of these systems was an increased 
interest in the overall effects resulting from the characteristics of the trans- 
mission path. This report deals with several of the losses associated with the 
transmission path including absorption, scattering, and turbulence. 

The bulk of the work has appeared earlier in the form of two NASA 
Technical Memoranda, "Atmospheric Transmission of CO 2 Laser Radiation 
With Application to Laser Doppler Systems" [7], and "Laser Doppler Systems 
in Atmospheric Turbulence" 1.8] . This work has been combined with an examina- 
tion of the propagation effects of aerosols in an attempt to more clearly define 
carbon dioxide transmission losses in an atmospheric path. 

Section n of this report discusses transmission losses by addressing 
those losses attributed to absorption, particularly absorption by carbon dioxide 
water vapor, and nitrous oxide. In addition, absorption losses are calculated 
for a 1973 flight test of the system at Edwards Air Force Base. A brief 
comparison of hydrogen fluoride (HF) and deuterium fluoride (DF) lasers with 
CO 2 lasers is also provided. Section in is a discussion of turbulence effects, 
including turbulence parameters and optical effects, and their effects on laser 
beam propagation and coherent heterodyne systems. It concludes with the cal- 
culation of the signal-to-noise equation. Section IV is involved with aerosol 
effects on the backscatter and extinction coefficients and their relation to the 
signal-to-noise ratio. The performance of the MSFC Pulsed Laser System at 
ground level is also discussed. The conclusions of this study are presented in 
Section V. 


II. ABSORPTION LOSSES 


Transmission losses are characterized as linear or nonlinear to dis- 
tinguish between laser power independent and dependent processes. In this 
section we examine the linear, power independent losses attributed to molecular 
gases and aerosols. 

An atmospheric window is defined as a spectral region over which the 
absorption of a particular wavelength by the permanent atmospheric gases is 
relatively weak. Although none of the windows are completely transparent 
because of molecular absorption and particulate scattering, selection of a 
window is mandatory for long-range propagation. The atmospheric windows in 
the infrared wavelengths lie between strong bands of carbon dioxide and water 
vapor, the two major absorbing gases of infrared. Figure 1 shows a low resolu- 
tion absorption spectrum (1 to 15 /Ltm) of solar radiation at ground level [9]. 

In addition, the spectra of various atmospheric gases are displayed as obtained 
from laboratory measurements. The operating wavelength regions of HF, DF, 
and CO 2 lasers are indicated in the figure. The gases indicated (with the excep- 
tion of ozone and water vapor) are generally uniformly mixed. It should be 
noted that within the infrared region relatively narrow line widths can be 
obtained with CO 2 lasers. These lines are typically 2 x 10~® pm or 60 MHz, 
while the absorption lines of the atmospheric gases are several orders of mag- 
nitude wider. The transmission losses for several CO 2 lines are shown in 
Table 1. These values were obtained for a 10 km horizontal roundtrip path at 
an altitude of 3. 5 km. While the P40 and RO lines have the least transmission 
loss, practical considerations involving availability of equipment have necessi- 
tated the choice of 10.6 pm for operation of the MSFC system. McClatchey and 
and Selby [10] have calculated the atmospheric transmission in the 10.6 pm 
region for a set of standard atmospheric conditions. Their results are shown 
in Figure 2 as round trip losses for several P and R lines using the Air Force 
Cambridge Research Laboratories (AFCRL) Midlatitude Summer Hazy 
Atmosphere. 


A. Absorption Coefficient 

Since atmospheric conditions are variable and deviations from the stan- 
dard models can result in large errors, a more accurate indication of absorp- 
tion loss can be obtained using the prevailing atmospheric conditions for a given 
e:q)eriment. This can be accomplished by providing an absorption coefficient 
for each of the molecules relevant to CO 2 laser radiation as a function of pressure 


2 


ABSORPTION 



Figure 1. Near-infrared solar spectrum, spectra of various 
atmospheric gases and operating frequencies of HF, DF, 

and CO 9 lasers. 





TABLE 1. TRANSMISSION LOSSES 


Laser Line 

Wavelength (pm) 

Loss (dB) 

P40 

10.81 

2.15 

P32 

10.72 

5.00 

P20 

10.59 

7.85 

P18 

10.57 

10.40 

P14 

10.53 

9.23 

P6 

10.46 

4.64 

P4 

10.44 

3.58 

RO 

10.40 

2.20 

R12 

10.30 

6.80 

R22 

10.23 

6.30 

R30 

10. 18 

4.60 


and temperature. The absorption coefficient as a function of frequency is 
assumed to be described by the Lorentz relation 


K(i^) 


So 

Tr[(u - U )^ + Q!^] 
o' 


( 1 ) 


where 


- the resonant frequency (cm 

S = the line intensity per absorbing molecule (cm“^ /molecules cm“^) 
a = the Lorentz line width parameter (cm“^). 


The validity of equation ( 1) in describing the correct line shape in the wings is 
in doubt, especially for carbon dioxide and water. In the first case, the wings 
are overpredicted, while in the latter, they are underestimated. 
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TWO WAY TRANSMISSION LOSS (dB/km) 



The line intensity, S, in equation (1) is independent of pressure but not 
temperature; i.e.. 


S(T) = S(Tg) 


Q (T ) 
V s' 


Q (T ) 
s' 


exp 


1.439 E”(T - T^) 

t"t 

s 


( 2 ) 


where 

T = the ambient temperature 

T = the standard temperature (296 K) 
s 


= the vibrational partition function 


= the rotational partition function 
E” = the energy of the lower state of the transition (cni~^) . 


Returning again to equation (l), the half-width of the Lorentz line, a , 
is dependent upon both temperature and pressure; i.e, , 


P 

a - a ~ 
o P 

o 
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( 3 ) 


where P and T are the ambient pressure and temperature, and P^, and 

are reference values. The parameter n is usually assumed to be equal to 

1/2. Although Benedict and Kaplan [11] suggest n= 0.62 for water vapor, and 
measurements by Ely and McCubbin [12] and Tubbs and Williams [13] suggest 
a value closer to n = 1 for carbon dioxide, the exact dependence is uncertain 
and appears to rely on the rotational quantum number of the line. As a result, 
the more generally accepted value of n= 1/2 will be assumed in the calculations 
of this report. The remaining parameters encountered in the calculation of the 
absorption coefficient were obtained from McClatchey et al. [ 14] who have 
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listed the four parameters u^, S, a and E” for seven molecules-, water, H2O; 

carbon dioxide, CO2; Ozone, O3; nitrous oxide, N2O; carbon monoxide, CO; 
methane, CH4; and oxygen, O2. 

The particular molecules which are active in the 10 jum region of the 
infrared are: methane, CH4; ethylene, C2H4; nitrous oxide, N2O; ammonia, 

NH3; ozone, O3; nitric acid, HNO3; carbon dioxide, CO2; and water, H2O. Among 
these molecules, carbon dioxide and water vapor are the major atmospheric 
attenuators at the CO2 laser line frequencies. The primary reason for this is 
that the atmospheric concentrations of the trace gases and pollutants are too low 
to have any significant effect. This is illustrated by the following simple 
estimate: At the center of the absorption line, the absorption coefficient is 
given by K = S/na per molecule. For a concentration of 1 ppm, the number of 
molecules in a 1 km path is approximately 10^®. For S = 10“^^ an tto = 0. 1, one 
obtains K = 10“^ c/km, where c is the number of ppm. Most trace gases have a 
concentration, c, of unity or less. For example, ozone has a concentration of 
0.02 ppm at sea level and 0.2 ppm at 25 km. Nitrous oxide, in turn, has a con- 
centration of 0.28 ppm. Consequently, it is apparent that these effects will be 
very small compared to those of carbon dioxide and water. It should be noted, 
however, that at altitudes above 12 km, ozone absorption may become more 
important (particularly at the R-lines). Because of the small effect of the 
trace gases, calculations in this report are primarily restricted to carbon 
dioxide, water vapor, and nitrous oxide. The calculations for nitrous oxide are 
performed to indicate the level of trace gas contributions to the overall absorp- 
tion. It should be remembered that in addition to molecular absorption, scatter- 
ing and absorption due to aerosols must be considered. 


B. Absorption by Carbon Dioxide 

The absorption coefficient at the center of P20 laser line by the atmos- 
pheric carbon dioxide has been calculated by Yin and Ling [15] who have pro- 
vided polynomial fits to the absorption coefficient as a function of altitude for 
two model atmospheres. The wings of neighboring lines also contribute to the 
absorption at any laser line frequency. In this work, lines within ±20 cm“^ are 
assumed to contribute to the absorption. A concentration of 330 ppm of carbon 
dioxide in the atmosphere is assumed. The effect of including the wings of the 
neighboring lines is approximately an increase of 3 to 5 percent at the P20 line. 
Since the CO2 laser operates efficiently around the P20 line, the absorption 
coefficients for P16, P18, P20, P22, and P24 lines are calculated for tempera- 
tures and pressures of interest in the lower atmosphere. The effect of tempera- 
ture on the absorption coefficient may be seen from Figures 3 and 4. The 
absorption coeffieients for P16, P18, and P20 lines are approximately the same 
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ABSORPTION K (km'') 



Figure 3. Absorption coefficient of carbon dioxide at T = 300 K. 
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ABSORPTIOM COEFFICIENT K (km’’ ) 



Figure 4. Absorption coefficient of carbon dioxide at T 



at 300 K and they separate out at 200 K. The absorption coefficients are almost 
independent of pressure and the slight pressure dependence comes from the line 
wings as the absorption at the line center is independent of pressure. Figures 5 
through 9 give the absorption coefficient at the individual laser lines for pres- 
sures from 100 to 1100 mb and temperatures from 200 to 300 K. 


C. Absorption by Water Vapor 

Many weak absorption lines of water vapor occur in the 10 nm region. 

On the basis of these lines, there should be an almost complete transparency in 
this region for reasonable atmospheric water content. However, measurements 
of heat radiation of the sky yield higher values of total radiation than expected 
due to the rotational lines of water vapor. To explain the observed higher 
emissions, Elsasser [16] in 1938 postulated the existence of a water vapor con- 
tinuum in the 8 to 14 pm atmospheric window and beyond. This continuum was 
believed to be due to the wings of strong lines located in the bands on either side 
of the window. The existence of this continuum has now been well established 
as a result of several experiments in solar spectrum observation and laboratory 
measurements [17-22] . However, the theory based on usual line shapes has 
still not been successful in explaining experimental observations; consequently, 
calculations have to be based on available experimental results. 

McCoy, Rensch, and Long [22] measured the absorption of the water 
vapor equation and suggested the empirical equation 

K= 4.32 X 10"® p (P + 193 p) km"^ , (4) 


where p is the partial pressure of water vapor and P is the total pressure, 
both in torr. This equation solves the problem of pressure dependence and is 
valid at room temperature. This relation has been widely used to estimate 
water vapor absorption. 

Burch [23] measured the continuum absorption coefficient for pure water 
vapor at three temperatures and his results are shown in Figure 10. The points 
marked are revised values and are 10 to 15 percent less than the solid curves 
obtained in earlier measurements [24] . The rapid decrease in the absorption 
coefficient with increasing temperature is a prominent feature of these results. 
This trend was predicted by Varanasi, Chow, and Penner [25] as being a result 
of the association of water molecules due to hydrogen bonding (dimerization). 
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Figure 7. Absorption coefficient of carbon dioxide 
at V = 944. 190 cm“^ . 
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Figure 10. Continuum absorption coe 
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For atmospheric transmission studies, we are interested in temperatures 
less than 296 K. Since no data are available at low temperatures, an empirical 
extrapolation suggested by Burch is used in this work: 


(T) 


(296) 



( 5 ) 


where C^° is the absorption coefficient (molecule"^ cm^ atm~^) and t is to be 

determined from the experimental values. At 10.6 pm, we find that t= 5.300 
for T = 358 K and t = 5. 25 for T = 392 K. t appears to decrease slightly with 
temperature. These t values are for temperatures above room temperature; 
no laboratory measurements are available for temperatures below room tem- 
perature. One value of the absorption coefficient at T = 283 K is 10 gm~^ cm^ 
atm“^ from Figure 8 of Bignell's [21J paper from the solar spectrum observa- 
tions. The magnitude of the error is unknown here. This value corresponds to 
T = 5.1 at T = 283 K. Thus, the value of t appears to remain approximately 
5 even below T = 296 K. We assume t= 5.25 in this work. 

There is also the effect of foreign gas broadening on the continuum 

absorption. The results of McCoy, Rensch, and Long [22] give C ® = 0.005 C " 

B s 

at room temperature. The continuum absorption coefficient is given by K = Cu 
where u is the absorber thickness expressed in molecules cm“^ and C is given by 


C = C °p + C ®P 
s ^ B 


( 6 ) 


where p is the partial pressure of water vapor and P is the total pressure, 
both in atmospheres. The partial pressure of water vapor may be obtained from 

p = 4.56 X 10"^ w T (7) 


where w is the number of gm-cm“^/km of water vapor. The water vapor is 
usually expressed in units of precipitable centimeters (pr-cm) which is the same 
as gm-cm~^. The absorber thickness u is given by 
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u 


7. 18 X 10^6 p , , . 9 

— molecules/cm^ 


( 8 ) 


= 3.34 X 10^2 


w = 


The continuum absorption may be written as follows after combining equation ( 5) 
with equation ( 8) : 


K 


1.58X 10^ / 296 ^'^^ 
T \ T / 


p(p + 0.005 P) km-‘ 


( 9 ) 


We used C ° = 2.2 x 10“^^ molecules ^ cm^ atm ^ at 10.6 pm. Equation (9) is 
s 

used in this work to calculate the water vapor continuum absorption coefficient 
for all CO 2 laser lines shown in Figure 11. 

Recently McClatchy and D'Agali [26] gave the following ejqjression for 

the temperature variation of C (t'jT) in the 8 to 14 pm continuum.- 

s 




C ( p, 296) exp 
s 


1800 


(7 



( 10 ) 


where 


C ( 1^,296) = 4. 18 -f 5578 exp (-7. 87 X 10"3 i^) 
s 


( 11 ) 


and Cj^(v,T) = 0.002 C^(i/,T) . The unit of T) and C^( t', T) is (pr-cm) “1 

atm“^. For temperatures within ±15° of 296 K, the power law variation in equa- 
tion (5) and the exponential variation in equation (lO) yield absorption coefficients 
within 11 percent of each other. However, the exponential dependence has better 
theoretical justification and equation (10) is preferable for the water vapor 
continuum. 

In addition to the continuum, rotational lines of water vapor absorb in 
the 10.6 pm region. Line-by-line calculation is performed for lines which are 
within ±20 cm~^ from the laser line. The line parameters tabulated by McClatchy 
et al. are used in these calculations. The effects of varying the pressure and 
the temperature are shown in Figures 12 through 18. 
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D. Absorption by Nitrous Oxide 

Nitrous oxide has several lines in the 10 region. The absorption due 
to these lines at the laser line frequencies is obtained by a line-to-line calcula- 
tion assuming that lines within ±20 cm“^ contribute. The concentration of nitrous 
oxide is assumed to be 0.28 ppm. The results are shown in Figures 19 to 25 for 
various pressures and temperatures. It may be seen that the absorption of nitrous 
oxide is approximately two orders of magnitude less than that of carbon dioxide. 


E. Comparison with AFCRL Calculations 

The present calculations are compared in Table 2 with those of McClatchey 
and Selby [10] for the AFCRL Standard Atmospheric Models at sea level. Our 
results have been within ±20 percent of AFCRL calculations. 


F. Absorption Losses for Edwards AFB 1973 
Flight Test 

The 1973 flight test of the MSFC Pulsed Laser Doppler System occurred 
aboard the NASA CU990, Galileo I. The particular test in consideration involved 
the descent of the airplane against the uniform target of Roger’s Dry Lake on 
Edwards AFB. During this test, the atmospheric temperature, dew point, and 
altitude were measured along the flight path in addition to the particle concentra- 
tion and laser system signal -to -noise (S/N) ratio. The data from this test are 
used to calculate the transmission loss for five P lines and assess the effect on 
the S/N. A detailed derivation of the S/N, including the effects turbulence, is 
discussed in a later section. In this case, the simplified expression for the S/N 
is given by 


S/N 


E a’ 

I6hv[L»+(^) 


( 12 ) 


The parameters involved for this calculation are as follows [27]: 
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Figure 22. Absorption coefficient of nitrous oxide 
at V = 942. 38 cm~^ . 
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Figure 23. Absorption coefficient of nitrous oxide 
lines at V = 944. 19 cm“^. 
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Figure 25. Absorption coefficient of nitrous ojdde 
lines at V = 947.738 cm“^. 






TABLE 2. COMPARISON WITH THE RESULTS OF 
McCLATCHEY AND SELBY 



Tropical 

Midlatitude 

Summer 

Midlatitude 

Winter 

Subarctic 

Summer 


AFCRL 

Present 

AFCRL 

Present 

AFCRL 

Present 

AFCRL 

Present 

P16 

0.5722 

0.5481 

0.3527 

0.3728 

0.07466 

0.08775 

0.1953 

0.2274 

P18 

0.6346 

0.4972 

0.4146 

0.3377 

0.1223 

0.0907 

0.2548 

0.2065 

P20 

0.6094 

0.5051 

0.3852 

0.3419 

0.09575 

0.08964 

0.2238 

0.2078 

P22 

0.6058 

0.4842 

0.3867 

0.3266 

0.1021 

0.084 

0.229 

0.1973 

P24 

0.6029 

0.4796 

0.3815 

0.3216 

0.09554 

0.0795 

0.2223 

0. 1921 


Parameter 

Signal -to-Noise Power Ratio 

Laser Pulse Energy 

Modified Target Cross-section 

Target Range 

Aperture Diameter 

System Efficiency 

Detector Quantum Efficiency 

Atmospheric Transmission 

Wavelength 

Energy per Photon 


Symbol 

Value 

S/N 

E 

12 mj 

CT* 

L 

0.05 sr“^ 

d 

10 in. 


0.05 



1.06 X 10 

hv 

1.9 X 10"' 


The atmospheric transmission loss is calculated from 


V 


a 


L 

exp / 2 K (L) dL 
0 


(13) 


The range is converted into altitude by the relation 


H = 1.97 kft+ L sin 7“ 


(14) 
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The loss in S/N ratio in dB due to atmospheric absorption at any range is 
calculated from 



atm loss 


H 

8.68 f K(H) 
1.97 


dH 

sin 7° 


(15) 


The measured values of the static temperature (T^) and the dew point 

temperature (T^^) along the flight path at various altitudes for flight B8, Run 18, 

1/19/1973, at Edwards AFB are given in Table 3. The individual molecular 
absorption coefficients corresponding to these atmospheric conditions are calcu- 
lated from the data of this report and are shown in detail for the P20 line. The 
aerosol attenuation must also be added to obtain the total absorption coefficient. 
The aerosol absorption and scattering coefficients for the hazy condition are 
taken from Reference 10 and are listed in Table 4 together with the total molecu- 
lar absorption coefficients for five P lines. The atmospheric two-way trans- 
mission loss is calculated using equation (15) and is shown in Figure 26 for the 
five P lines. The P22 and P24 lines have lower loss and P16 and P18 lines have 
higher loss than the P20 line. The magnitude of the difference is approximately 
2 dB over a path of 90 kft. Thus, the variation from line to line is not significant 
around the P20 line. The S/N ratio corrected for the transmission loss is shown 
in Figure 27 together with the measured values. If there were no atmospheric 
losses, the curve falls off as 1/R^ where R is the range. The measurements 
indicate a fall -off nearly as l/R®. Though the CO 2 laser operates in the atmos- 
pheric window, the transmission loss is significant over long paths which may 
explain a large part of the losses leading to l/R^ dependence of the measured 
S/N values. Some of the remaining losses may be of instrument origin and 
some due to the atmospheric turbulence. 

The range of the system will improve if the large transmission loss can 
be reduced. One way to reduce the transmission loss is to operate at a different 
infrared wavelength having less atmospheric absorption. This possibility is 
examined in the next section. 


G. Lasers with Low Transmission Loss 

Infrared transmission of HF and DF laser radiation has been calculated 
by Wang [28], Spencer, Denault, and Takimoto [29] have given experimental 
results for several DF laser lines. We will now use this information to calcu- 
late the transmission loss for the laser Doppler system using HF and DF laser 
lines and compare with the CO 2 laser system. 
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TABLE 3. CALCULATION OF ABSORPTION COEFFICIENT FOR FLIGHT B8, 

RUN 18, 1/19/1973, EDWARDS AFB 


H (kft) 


T “C 
0 

P20 

CO 2 

P20 

H 2 O 

Continuum 

P20 

H 2 O 

Lines 

P20 

N 2 O , 

P20 

Total 

km"* 

2.39 

-0.6 

-7.9 

0.05187 

0.02773 

0.000983 

5 X 10“® 

0.08063 

3.00 

-3.6 

-9.4 

0.04821 

0.02633 

0.000818 

3.9 X 10"® 

0.0754 

4.00 

-8.6 

-9.9 

0.0424 

0.0309 

7.1 X 10"^ 

3.5 X 10"® 

0.07405 

5.03 

-10.5 

-11.1 

0.0404 

0.0291 

6.08X 10"^ 

3.3 X 10"® 

0.07014 

5.98 

-9.9 

-13.8 

0.04105 

0.02016 

4.75 X 10“^ 

2.9*x 10"® 

0.0617 

6.98 

-12.3 

-15.3 

0.03852 

0.01877 

3.95X l0-‘‘ 

2.6 X 10"® 

0.05771 

8. ob 

-13.5 

-16.6 

0.03729 

0.0168 

3.37X lO"'* 

2.4 X 10"® 

0.05445 

8.99 

-15.0 

-19.0 

0.03579 

0.01346 

2.61 X lO'-^ 

2.1 X 10"® 

0.04953 

10.00 

-16.4 

-20.8 

0.03444 

0.0115 

2.13 X 10-^ 

1.9x 10"® 

0.04617 

11.03 

-17.3 

-25.9 

0.03358 

0.0068 

1.29 X 10“'* 

1.8 X 10"® 

0.04053 

12.02 

-18.5 

-27.4 

0.03246 

0.00568 

1.07 X 10"'* 

1.6 X 10"® 

0.03826 

12.96 

-19.5 

-34.2 

0.03156 

0.00267 

5 X 10"® 

1.5 X 10"® 

0.0343 

14.00 

-19.8 

-45.0 

0.03128 

0.00075 

2 X 10"® 

1.3 X 10"® 

0.03205 

14.99 

-20.7 

-42.3 

0.03048 

0.001044 

2 X 10"® 

1.2 X 10“® 

0.03156 

16.01 

-21.3 

-40.6 

0.02997 

0.00127 

2.4 X 10"® 

9.5 X 10-® 

0.03127 


CO 



CO 

00 


TABLE 4. MOLECULAR ABSORPTION COEFFICIENT OF P LINES AND 
AEROSOL ATTENUATION PER KILOMETER 


H (kft) 

K (P16) 

K (P18) 

K (P20) 

K (P22) 

K (P24) 

K + cr 
a a 

(aerosol) 

2.39 

0.08803 

0.08171 

0.08063 

0.07512 

0.07062 

0.0248 

3.00 

0.08227 

0.07658 

0.0754 

0.07032 

0.06603 

0.02075 

4.00 

0.08081 

0.07542 

0.07405 

0.06963 

0.06572 

0.0153 

5.03 

0.07634 

0.07143 

0.07014 

0.06585 

0.06207 

0.0114 

5.98 

0.06717 

0.06303 

0.0617 

0.05747 

0.05363 

0.0085 

6.98 

0.06278 

0.05903 

0.05771 

0.05374 

0.05007 

0.0057 

8.00 

0.05894 

0.05578 

0.05445 

0.05062 

0.04703 

0.0036 

8.99 

0.05386 

0.0508 

0.04953 

0.04588 

0.0442 

0.0026 

10.00 

0.05024 

0.04751 

0.04617 

0.04267 

0.0393 

0.002 

11.03 

0.04378 

0.04157 

0.04053 

0.03686 

0.03352 

0.0017 

12.02 

0.04167 

0.0396 

0.03826 

0.035 

0.03175 

0.0013 

12.96 

0.03736 

0.03563 

0.0343 

0.03125 

0.02796 

0.001 

14.00 

0.03487 

0.03339 

0.03205 

0.029 

0.02579 

0.00072 

14.99 

0.03438 

0.03288 

0.03156 

0.02854 

0.02542 

0.00053 

16.01 

0.03409 

0.03258 

0.03127 

0.02829 

0.02521 

0.00039 








Figure 27. Comparison of measured and theoretical 
S/N values. 


The performances of CO 2 , HF, and DF laser Doppler systems looking 
down at a ground target from 5 km altitude at an inclination of 15 deg in AFCRL 
Midlatitude Summer Atmosphere are shown in Figure 28. The curves without 
atmospheric effects are calculated for the same values of the parameters 
assuming that the backscatter coefficient varies as lA^ [30] . The CO 2 laser 
experiences a loss which is 16 dB higher than that for DF laser for a slant 
range of 20 km. 
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Figure 28. Performance of CO 2 , HF, and DF laser Doppler systems against 
ground target at 15 deg inclination in AFCRL Midlatitude 
Summer Atmosphere. 
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The round-trip transmission losses for horizontal transmission over a 
range of 40 km at altitudes of 5 km and 10 km are given in Table 5. 

The DF laser has the least loss. As the altitude increases, the magni- 
tude of the loss reduces. At 10 km altitude, HF, DF, and CO lasers have 
approximately the same loss. While DF laser is much better than the CO 2 laser 
for transmission through the atmosphere, due to the longer wavelength the CO 2 
laser system has a larger coherence diameter which may be an important con- 
sideration for an operating system. 

TABLE 5. ROUND-TRIP HORIZONTAL TRANSMISSION LOSS FOR 
INFRARED LASERS OVER 40 km RANGE 



Loss in dB 

Laser, Line, and Wavelength 

H = 5 km 

H = 10 km 

CO 2 , P20, 10.6 pm 

19 

7 

HF, P 2 ( 8 ), 2.91pm 

4 

0.6 

DF, P 2 ( 8 ) , 3. 8 pm 

1 

0.45 

CO, Pg(5) , 5. 07 pm 

3.5 

0.45 
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III. TURBULENCE LOSSES 


The atmospheric turbulence resulting from random temperature fluctua- 
tions causes a variety of optical effects: twinkling (variation of image bright- 
ness) , quivering (displacement of image from normal position) , tremor disk 
( smearing of the diffraction image) , dancing ( continuous movement of a star 
image about a mean point) , wandering (slow oscillatory motions of the image 
with a period of approximately 1 minute and angular excursions of a few seconds 
of arc) , pulsation or breathing (fairly rapid change of size of the image) , boiling 
(time-varjdng nonuniform illumination in a larger spot image) , and image dis- 
tortion [31]. In the case of the laser beam, the turbulence effects are broadly 
categorized as beam wander, spot dancing, beam spread, and scintillations. 
These effects are caused by the temporal and spatial fluctuations of the direction, 
phase, and intensity of the wavefront. A description of the atmospheric 
inhomogeneities causing these effects is giveninthis section. 


A. Parameters of Atmospheric Turbulence 

The optical effects of interest are produced by the variations of refrac- 
tive index along the path of the beam. The changes in the refractive index are 
caused by the fluctuations of temperature which arise in turbulent mixing of 
various thermal layers. It is found that the temperature fluctuations obey the 
same spectral law as the velocity fluctuations. In analogy with the velocity 
turbulence, the atmosphere may be imagined to consist of a large number of 
eddies with varying dimensions and refractive indicies. For isotropic tur- 
bulence, the velocity spectrum 0^(K) (where the wavenumber K= 2ir/2 , i 

being the size of a turbulent eddy) contains the information about the turbulence. 
0^(K) is intuitively interpreted as the amount of energy in the turbulent eddies 

of size i . The kinetic energy of turbulence is assumed to be introduced through 
eddy scale sizes larger than the "outer scale" of turbulence, L^, corresponding 

to a spatial wavenumber = 2ir/'L^. As the wavenumber increases beyond K^, 

the turbulence tends to become isotropic and homogeneous. Experiments have 

11/3 

confirmed the predictions of Kolmogorov's theory of K dependence for the 

three-dimensional spectrum up to wavenumber K related to the inner scale of 
_ m 

turbulence by K =5. 92/i„. The spectrum in the dissipation region for 
0 m 0 

which the scale sizes are smaller than or wavenumbers greater than K 

0 m 
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is steeper than K . The spatial wavenumber region between the inner scale 
and the outer scale obeying the Kolmogorov theory of a spectral slope of 

-11/3 is known as the inertial subrange, is of the order of several milli- 
meters and Lq several meters. 

At optical wavelengths the refractive index may be related to the tem- 
perature throu^ the relation 


77.6 pf. 0.00753 

ni = n - 1 - — j^l + ^^2 X 10 J 

where 

n = refractive index 
P = pressure (mb) 

T = temperature (K) 

A = wavelength (/jm) . 


(16) 


The refractive index structure constant is a measure of the strength of 

atmospheric turbxilence and is defined as the mean-square difference in the 
refractive index at two points divided by the separation distance r raised to the 
2/3 power, i.e. , 


2 _ ^(ti9 - n^)^^ 

n ” 


(17) 


The temperature structure constant, , is similarly defined for tem- 
perature fluctuations. Thus, may be determined from temperature measure- 
ments alone from the following relation: 


C 

n 


77.6 P 

rji2 



0.00753' 

a2 


X 10“® c 

T 


(18) 
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Typical temperature variations of approximately 1 K resulting in refrac- 
tive index variations of a few ppm will cause significant effects on optical radia- 
tion fields propagating through the atmosphere. is measured in units of 
- 2 / ^ 

(meters) ^ and varies from 10”*® or more for strong turbulence to 10”** or less 

for weak turbulence. 

Assuming that the form of the temperature spectrum and the refractive 
index spectrum is the same as that of the velocity spectrum, Tatarski [32] gave 
the following form for the refractive index spectrum: 





( 19 ) 


where = 5.92. This spectrum has a singularity at K = 0. Physically, 

the singularity implies that the energy per unit volume becomes unbounded as 
the eddy size increases. This singularity is avoided by using the modified von 
Karman spectrum [33]: 





( 20 ) 


which gives a bounded variation for K < 2tt/Lq. The three-dimensional refrac- 
tive index spectrum models according to equations (19) and (20) are shown in 
Figure 29. Deviations from these models occur if the turbulence is not homo- 
geneous. Deviations have been observed in measurements over paths close to 
the ground and under conditions of weak turbulence [33] . 


B. Optical Properties of the Turbulent Eddies 

A qualitative description of the effect of the optical inhomogeneities on 
the wavefront is sometimes helpful. As a plane wavefront passes through a 
region of relatively lower refractive index, the wavefront will advance as the 
speed of light is higher. When passing through a region of higher refractive 
index, the wavefront is slowed down. Thus, the refractive inhomogeneities of 


45 



SPECTRUM AMPLITUDE 



R = 27T/C WHERE C- SIZE OF TURBULENCE EDDY; 
Kg = 27 T/Lq, typically Lq - 1 - 100 m; 

“ Ztt/Lq, typically Cq = 1 - 3 mm 
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Figure 29. Three-dimensional turbulence spectrum models 
for index-of-refraction fluctuation. 


the atmosphere deform the wavefronts. Portions of the wavefront maybe con- 
vex or concave in the direction of propagation. The concave portions converge 
while the convex portions diverge. Thus, it is possible for the intensity to 
build up in some areas. Since the inhomogeneities are random and moving, the 
bright and the dark areas move randomly. This point of view has been developed 
by the optical astronomers [31] . 

In general, the optics of the turbulent atmosphere may be described in 
teiTiis of a collection of weak, moving, three-dimensional, gaseous lenses of 
varying scale sizes bounded by the inner scale io the outer scale Lq of 
the turbulence. From equation (17) , it may be seen that the refractive index 
fluctuations associated with any scale size i increases as the size of the eddy 
increases: 


<[6n(i)]2> 



n 


( 21 ) 


An eddy of radius i can act as a lens with a focal length f = R/n^, where 
R is the radius of lens curvature. This is approximated from equation (21) as 


f 



i 



( 22 ) 


Turbulence as it naturally occurs can consist of eddies of hot and cold 
air. Eddies cooler than the ambient air have a higher index of refraction. In 
this case, nj is positive and the lens is convergent. If the eddy consists of 
warmer air, the index of refraction is lower (giving negative n^) and the lens 
is divergent. Typical values in the atmosphere are ~ 10“\ In^l ~ 10~®, 

and i a few centimeters. Thus, the focal length of the atmospheric lenses 
is f ~ 10 km. 

The turbulent eddies cause diffraction and refraction of the optical wave- 

front, and it is instructive to consider these two effects separately. Each eddy 

of size i bends as a ray bundle individually and has a diffraction pattern of 

approximate width, XL/S. , where L is the distance from the eddy to the receiver. 

The diffraction effects are negligible if XL/S « S or if i >> n/ XL,' and this is 

the geometric optics regime. But, if XL/S » S , the diffraction pattern entirely 

determines the intensity distribution of the receiver. Thus, the distance L . 

d 

beyond which the diffraction effects are important may be defined by ~ S^/X. 
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The refractive bending of rays by an eddy of size i depends on the 
refractive index of the eddy and may be calculated by Snell’s law to give 60 
6n. 6n is zero on the average. The mean square angular fluctuation for one 
eddy is <(50)^> ~ <(6n)^>. 

Over a path length L, let the beam traverse L/i eddies. Since the 
deflections are uncorrelated from one eddy to another, they add up as mean 
squares and the total mean square deflection is 

. ( 23 ) 


The effect of refraction is to change the cross-section of the beam from 
to (i + L0 )^. The amplitudes of the wave before and after refraction are 
related as 


= (A + AA)^ (i ± L0^^ 


(24) 


The fractional amplitude change is then given by 


X = 



(25) 


The variance of the amplitude fluctuation is then given by 




-V 

C^L^i /3 
n 


(26) 


The particular scale size 1 which contributes predominantly to the amplitude 
fluctuation depends on the path length; i.e. , different eddies are more effective 
at different path lengths in producing amplitude fluctuations. This point has 
been discussed by Tatarski [34] and De Wolf [35] as follows. 
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There are no diffraction effects in equation (26) . Diffraction effects are 

2 

negligible if the path length L is short or L « £ /X in the geometric optics 

S S V/ 

regime. Since (inner scale length) is the smallest eddy size, the diffraction 

by eddies of all sizes is negligible in this case. The smallest eddies then have 
the largest contribution. Equation ( 26) then becomes 



n s 0 


for 


L « 
s 



(27) 


2 , 

For path lengths greater than L^, such that » £q/'^> the eddies 
with sizes greater than £^ (£ > £ ^) have negligible diffraction, while eddies with 
sizes smaller than £ ^ {£ < -^^) diffract the rays where 


= ^/al7 . 

£ £ 


The diffraction pattern entirely determines the intensity distribution and the 
refraction is negligible for the small eddies of size £ « £^> i*®* » these eddies 

simply scatter the energy, not focusing the rays, and have little effect on the 
amplitude fluctuations. Geometric optics results may still be applied if the 
eddies of size £ < £ ^ are excluded. Then the smallest eddy size contributing 

significantly to the amplitude fluctuation is £ ^ and £ - £ ^ is substituted in 
equation (26) : 




C^L„^(AL J 
n £’ 


■Vfi 


C-'K 

n 


1 / 11/ 
/Gt ^6 


for 


» 


£IL 


(28) 


This is also the result obtained by the Rytov approximation. Equation (10) has 
been interpreted by Tatarski [9] as the result of the application of geometrical 
optics to a medium from which all eddies of scale sizes less than \/ XL' have been 
removed. 
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Young [36] obtained equation (28) on considering an aperture of Fresnel- 
zone size while calculating the normalized variance or the total modulation 
power obtained with a circular aperture using geometric optics. He observed 

that ray optics gave the correct parametric dependence because very little power 

-1/ 

is contributed by the range of spatial frequencies greater than (AL) ^ which are 
in the wave optics regime. Thus, weak scintillation may be thought of as aper- 
ture filtering with a Fresnel -zone size aperture which filters out spatial fre- 

- 1 / 

quencies higher than (XL) 

A related view developed by radio astronomers is the phase screen con- 
cept. Immediately after passing through a layer of inhomogeneities, only the 
phase is affected and the intensity remains constant. As the distorted wavefront 
propagates further, the relative phases will be changed due to different directions 
of propagation and in this way, the amplitude fluctuations develop. This concept 
has been applied by Lee and Harp [37] to wave propagation by dividing the three- 
dimensional refractive field of the medium into thin slabs perpendicular to the 
direction of propagation. The effects of the refractive inhomogeneities on the 
propagation of laser beams in the atmosphere is discussed in the following 
section. 


C. Effect of Atmospheric Turbulence on Laser Beams 

A qualitative description of the effects based on the relative sizes of the 
beam and the refractive inhomogeneity is given first, followed by the results of 
Rytov approximation on the wave equation. 

If the scale size of the inhomogeneity is much larger than the diameter 
of the laser beam, the entire beam is bent away from the line-of-sight and 
results in beam wander or beam steering. The center of the beam executes a 
two-dimensional random walk in the receiver plane and, for isotropic turbulence, 
the displacement of the beam from the line-of-sight will be Raleigh-distributed 
[38] . Inhomogeneities of the size of the beam diameter act as weak lenses on 
the whole or parts of the beam with a small amount of steering and spreading. 
When the inhomogeneities are much smaller than the beam diameter, small 
portions of the beam are independently diffracted and refracted and the phase 
front becomes corrugated. The propagation of the distorted phase front causes 
constructive interference over some parts of the receiver and destructive inter- 
ference over others, leading to alternate bright and dark areas. Since the 
atmosphere is seldom stationary, the locations of the bright and dark regions 
change continually and give the pattern of "boiling.” Thus, the turbulence causes 
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the received field to scintillate in time and space. Since the atmosphere con- 
sists of inhomogeneities of all sizes in motion, the laser beam experiences 
fluctuations of beam size (spreading) , beam position (steering or wander) , and 
intensity distribution within the beam (scintillation) simultaneously. The rela- 
tive importance of these effects depends on the path length, strength of turbu- 
lence, and the wavelength of the laser radiation. The variance of irradiance 
reflects the deviations from the mean due to beam wander, broadening, and 
scintillations. Estimates of these effects for finite beams are discussed as 
follows. 

Finite transmitter effects have been treated by the Huygens -Fresnel 
approach [39] and by the transport approximation, Fante [40]. The Huygens- 
Fresnel principle was extended by Lutomirski and Yura [39] so that the second- 
ary wavefront is determined by the envelope of spherical wavelets from the 
primary wavefront as in vacuum, but each wavelet is determined by the propaga- 
tion of a spherical wave in the refractive medium for small scattering angles. 
Hence, a knowledge of the spherical wave propagation in the turbulent medium 
is sufficient to determine the field distribution at an observation point. The 
mean irradiance distribution at the observation point is characterized by the 
mutual coherence function (MCF) of a spherical wave and the geometry of the 
transmitting aperture. The MCF, defined as the cross-correlation function of 
the complex fields in a direction transverse to the direction of propagation, 
describes the loss of coherence of an initially coherent wave propagating in a 
random medium. When the medium is characterized by a refractive index 
variation that is gaussian with zero mean, the MCF is given by [32] 


M(p,L) = exp 




(29) 


where D is the wave structure function, p is the lateral separation between the 
two points in a plane perpendicular to the direction of propagation, and L is the 
path length of propagation. For spherical wave propagation and a Kolmogorov 
spectrum of turbulence [41] , the wave structure function is given by 


D(P) 


5 / ■*-' 

2.91 kV'^ / 
0 


ds 


C^(s) 
n ' 



k = 27 t/A 


(30) 
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For constant C 

n 


D(p) 


1.089 k^C ^Lp 
n 


% 


( 31 ) 


Equation (29) may be written as 


M(p,L) = exp 


■("f 


where 


Po = 


, io « P « L| 


0 » 


1.455 f ds C ^(s) 
0 


(r) 




(32) 


(33) 


For constant C 

n 


- 3 / 

Po = [0.545 k^LC2] ^5 


n 


(34) 


Po is the lateral separation such that the MCF becomes equal to l/e and is 
called the lateral coherence length. If d > po where d is the diameter of the 
aperture, the turbulence along the path of the beam reduces the lateral coherence 
between different elements of the aperture and effectively transforms it into a 
partially coherent radiator of dimension, pq, which decreases with increasing 
distance from the aperture [38] . But if d < po, the entire aperture behaves like 
a coherent radiator. Thus, the resulting field distribution is characteristic of a 
coherent aperture with a diameter d^^^ equal to the smaller of po and d. Yura 

[42] showed that if L >> kd d the beam properties may be approximated by 

elf 

those of a spherical wave and if L « kd beam properties may be 

approximated by those of a plane wave. 

In the absence of turbulence, a gaussian laser beam will have a far-field 
angular spread B - 2A./7T'd. When turbulence is present along the path, the 
scattering of the beam by the moving turbulent eddies causes additional spread- 
ing which is greater than 9 under some conditions. 
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The spreading of the beam consists of the deflection of the beam as a 
whole by the eddies larger than the beam diameter d and by the broadening of 
the beam due eddies smaller than d. The separation of the beam spread into 
the two components of wander and broadening is usually based on the length of 
the time of observation. If V is the transverse wind-velocity component, the 
turbulence eddies are interacting with the beam in time intervals of order 
At ~ d/V. For times much longer than At, the received spot will be fully 
spread due to wander and broadening. The total spot size, viewed on a time 
scale T » At, is called the long-term spot size and the beam radius of long- 
term spread is denoted by p . If a picture with an exposure time much less 

L 

than At were taken of the beam, a laser spot of radius p^ broadened by the small 
eddies at some distance p^ normal to the line-of-sight would be observed, p^ is 
the radius of beam-centroidal motion. 

The concepts of the beam spread just discussed apply only if the beam is 
intact. When the beam wander effect is removed, the atmospheric spreading of 
the beam at the receiver is the short-term beam spread; however, when the 
turbulence along the path is strong, the beam breaks up into multiple patches 
and there is not much wander. Kerr and Dunphy [43] found experimentally at 
near field that the focused beam broke up at the receiver into a proliferation of 
transmitter-diffraction-scale patches when there is either transmitter mis- 
adjustment or strong turbulence. Reidt and Hohn [44] observed in recent 
experiments that the focused beam at 0.63 pm was broken into several spots of 
characteristic length 2AL/7rd which is the free-space diffraction-limited focused 
radius. They verified also the inverse variation of the patch size with the 
diameter of the beam, that is, the increase in patch size with decreasing beam 
diameter or vice versa. 

The beam will be quivering while remaining intact as long as the turbu- 
lence is relatively weak and the path is not too long or, more precisely, if 

L < L where is the critical range defined by L = k d^ . where d is 
CR CR CR 6II 611 

the smaller of pg and d [20] . 

Quantitatively, the spreading of the beam depends on the lateral coherence 
length Pg. If the beam diameter d is much smaller than pg, the effect of turbu- 
lence on the beam is negligible and the beam diameter in the far field will be 
that due to diffraction, namely, 2AL/7rd. However, if the initial beam diameter 
is greater than the lateral coherence length pg, the beam diameter in the far 
field is given by 2AL/7rpg [45]. That is, the beam behavior is characteristic of 
an aperture of diameter pg instead of d. The turbulence along the path transforms 
the aperture into a partially coherent radiator. For the initial field distribution 
of gaussian form given by 
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, . r 2r^ ikr^ 

Uo(r) = exp — 


( 35 ) 


where F is the radius of the curvature, the long-term beam radius is obtained 
as [46] 



4L^ ^ 4JJ 

k^d^ ^4 \ F/ ^ 


(36) 


When the turbulence is strong, the beam will be broken up and for this case 
Klyatskin and Kon [47] found that the wander of the center-of-gravity of the 
beam may be obtained from 






(37) 


for constant C^, L » kp^ (far-field case), and d » p^. Taking the ratio of 

equations (37) and (3ft) and observing that the last term on the left-hand side of 
equation (36) is dominknt for strong turbulence, we obtain 



(38) 


Since L » kp^ then 


« 


which means that the motion of the beam 


centroid is negligible compared with the beam spread when the beam breaks up 
into multiple patches. It is not possible to tell how many patches will form, but 


the bright patches will be in a zone with mean square radius p [48], 

JLi 


For ranges less than the critical range, tracking laser transmitters have 
been proposed and discussed [49, 50] to compensate the beam wander which is a 
relatively low-frequency phenomenon of the order of 0. 5 Hz. For L < the 

beam centroidal motion is given approximately by [51] 
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( 39 ) 



2.97 iJ 


and appears to satisfy a gaussian distribution [52] on each axis. 

One of the basic observations on laser beams is the fluctuations in inten- 
sity in addition to fluctuations in beam size and position. Fluctuations in inten- 
sity are called scintillations. Experimental data on scintillations are extensive 
at present and are usually expressed as the normalized variance of the intensity; 


<I^> - <L>2 


(40) 


I is the measured value of the fluctuating intensity and is the average. The 
statistical theories of optical propagation in a randomly inhomogeneous medium 
usually calculate the variance in the log-amplitude X. The log-amplitude X(x) 
at any location x is related to the intensity or the irradiance l(x) by 


X(x) 




l(x) 

<o 


(41) 


where A is the amplitude. The variance of the log-amplitude is theoretically 
given by 


7 / 11 / 

oL = 0.124 
T n 


L » 


ii 


(42) 


for a spherical wave. Experiments have confirmed the increase of the log- 
amplitude variance as the 11/6 power of the path length, as the 7/6 power of 
the optical wavenumber k, and the square of the strength of turbulence for a 
homogeneous path, i.e. , = constant along the path as long as cr^ < 0. 3 for 

the spherical wave [53]. 
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Increase of the path length or the strength of turbulence beyond =0.3 
fails to produce corresponding increase in the measured value of the log- 

n 

amplitude variance a which slows down, reaches a maximum, and begins to 

decrease, reaching an asymptotic value of approximately 0.23 independent of 
wavelength, path length, and parameters of the turbulent medium. Thus the 
scintillation of an optical beam does not increase indefinitely but tends to 
saturate for strong turbulence and long paths. This is shown schematically in 
Figure 30. The more precise way of showing the saturation effect is to plot the 
theoretical value of log-amplitude variance along the x-axis and the experimental 
value of the log-amplitude variance on the y-axis, as is usually done in the 
literature [53]. 

The probability distribution of the amplitude or intensity fluctuation is an 
important aspect in optical communications. Theoretically the application of 
the central limit theorem leads to a normal distribution of the log-amplitude; 
i.e. , the scintillations follow a log-normal distribution. For a unit-amplitude 
plane wave, the probability density of intensity p(l) satisfies 


p(I) = n exp -(inl + ^j (2a2)-i 

(27r)/2 al L V 2/ 

where 

= In ^1 + j . 


(43) 


(44) 


For « 4, there is extensive experimental evidence in support of log-normal 

distribution of intensity. Physically, the reason for the log-normal distribution 
is the multiplicative effect of the refractive -index variations along the path on 
the intensity [34] ; i.e. , the field is modulated in each subrange of the path 
multiplicatively if it is assumed that the effect of atmosphere in each subrange 
is independent of the degree of coherence. Hence, the variation of the logarithm 
of the amplitude and intensity is the sum of several random variations induced 
along the path of propagation which leads to a normal distribution due to the 
central -limit theorem. Ej^erimental measurements by Gracheva, et al. [54] 
indicate that the log-normal distribution is a reasonable approximation for 
25 < cr^ < 100. For 1 < < 25, there are significant deviations from the 

log-normal distribution. They did not notice Rayleigh distribution, as predicted 
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■CALCULATED BY RYTOV METHOD 


■SATURATION REGIME 


EXPERIMENTAL DATA 


(PATH LENGTH) 


Figure 30. Saturation of the variance of log-amplitude fluctuations of a spherical wave 
propagating through a homogeneous ( = const. ) turbulence path. 




by de Wolf [55] for large values of Physically, Rayleigh distribution 

results when the field at the receiver is the sum of a large number of independ- 
ently scattered fields. The physical model of de Wolf [35] attributes the 
Rayleigh distribution to the multiple scattering by off-axis eddies which are 
large in strong turbulence. 

In addition to temporal variation, the received field exhibits a spatially 
varying intensity after propagation through turbulence. The spatial structure 
of scintillations is usually studied by the covariance of the logarithm of the 
amplitude fluctuations. A correlation length may be defined in several ways. 

One definition is the value of separation for which the spatial covariance function 
goes to zero for the first time. The correlation length provides a measure of 
the separation between two points over which the scintillations are correlated. 
The concept of the correlation length has been useful for some applications. 

For example, the theoretical calculation of the log-amplitude variance is 
usually based on a point detector. In practice, a collector aperture of diameter 
less than the correlation length is treated as a point detector for comparison of 
theory and experiment. Another use of the correlation length is in understanding 
the aperture averaging effect which is discussed later. The intensity covariance 
function is usually measured and is related to the log-amplitude covariance 
as [56] 


Cj(p) = <I> [exp (4 C^(p)} - 1] 


(45) 


where the separation is p = | x - x’ i > the covariance of log-amplitude is 


C^(p) = <[X(x) -<X>] [X(x') -<X>]> 


(46) 


and the covariance of intensity is 


Cj(p) = <[l(x) -<I>] [I(x’) -<!>]> 


(47) 


The intensity and the log-amplitude correlation distances may be seen to 
be the same from equation (45) but the curves differ greatly. For plane or 
spherical wave propagation, C (p) has the form 
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p 


(48) 


C^(P) - 




where f is a function which depends on the mode of propagation. The first zero 

1 / 

of f is the correlation length and is approximately 0.76 (A.L) ^ for a plane wave 
1 / 

and 1. 8 (XL) ^ for a spherical wave based on the perturbation theory of Tatarski. 

1 / 

Experiments confirm a correlation length of the order of (XL) ^ when there is 
no saturation scintillations. 


The classical method of reducing the fluctuations in the received power 
has been to increase the size of the receiving optics. The reduction of the fluc- 
tuations in the received power, when the diameter of the receiving aperture is 
increased, is called aperture averaging. The physical reasoning is that large 
apertures tend to average over the statistically independent portions of the 
scintillation pattern [56] . The diameter of an independent patch is usually 
imagined to be of the order of the correlation length. Thus for weak scintilla- 

tions, the size of the patch is of the order of the Fresnel length (XL) Experi- 
ments [57] show that the correlation length decreases as the turbulence gets 
stronger. Moreover, the covariance curve develops a progressively higher 
correlation tail and aperture averaging deteriorates as the turbulence strength 
increases. Figure 31 illustrates the covariance function for several 

To summarize, the phenomenon of saturation of optical scintillation may 
be characterized by the following experimental observations: 


1) There is a physical breakup of the beam into multiple patches when 
the turbulence is strong. 


2) The variance of the log-amplitude reaches a maximum instead of 
increasing indefinitely and then decreases to an asymptotic value of approxi- 
mately 0.23 independent of the wavelength of the optical radiation and parame- 
ters of the turbulent medium as the path length, or the strength of turbulence, 
or both are increased. 


3) The spatial covariance function of log-amplitude decays progressively 
faster but maintains a progressively higher tail as a function of the detector 
spacing with increasing turbulence strength and path length. 
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Figure 31. Theoretical curves of the log-amplitude covariance 
function in weak to strong turbulence [35] . 

4) Aperture averaging becomes progressively less effective for a given 
size of the aperture as the turbulence strength and the path length are increased. 

There were several theoretical attempts to explain the saturation 
phenomenon in the past but it is only in recent years that the theory is able to 
account for all the observed effects. Since the saturation effect is observed 
under conditions of strong turbulence, a working definition is needed for the 
conditions specified as weak or strong turbulence. The turbulence is called 
"weak” if the theoretical variance of log-amplitude 

= 0.124 k^®L (49) 

T n ' ' 
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satisfies the inequality < 0.3. The turbulence is called "strong” if cr^ > 0.3. 

Thus the specification of the turbulence conditions involves not just the refractive - 

index structure constant but a combination of the wavelength of radiation, path 

length, and . 

n 

One of the earlier attempts to physically explain saturation is that by 
Young [58] . When stellar scintillation is observed with small apertures, the 
total modulation power first increases with zenith angle and then decreases. 

Young suggested that the smearing or blurring of the details of the wavefront 
distortions caused by the intervening atmospheric turbulence is responsible for 
the observed saturation of astronomical scintillation. The blurring effect is 
produced by the small-scale fluctuations in the wavefront which change rapidly 
due to the short lifetime of the small eddies. Thus, a shadow pattern produced 
by a high layer in the atmosphere becomes smeared or washed out by the inter- 
vening turbulence as it propagates down through the atmosphere. 

Clifford, Ochs, and Lawrence [59] give a detailed description of the 
smearing process in the scintillation pattern in strong turbulence. The wave- 
front develops small-scale and large-scale distortions while passing through the 
turbulent medium. When such a distorted wave is incident on a Fresnel-zone- 
size eddy, the small-scale distortions reduce the resolving power of the eddy 
while distortions larger than the eddy merely tilt the diffraction pattern of the 
eddy. The small scale distortions may be expected to change in detail several 
times during the lifetime of the larger eddies due to the shorter lifetimes of 
small eddies. These fluctuations of the small-scale distortions of the wavefront 
cause the smearing of the diffraction pattern of the Fresnel-size eddy, thus 
making it less effective as a producer of scintillation. Scintillations remain 
bounded and do not build up indefinitely due to the smearing by the small-scale 
eddies. 


The smearing effect of strong turbulence is incorporated into the first- 
order analysis by means of a spectral -filter function which filters out higher 
spatial frequencies. The spectral filter function is given by the normalized 
two-dimensional Fourier transform of the short-term average of the irradiance 
profile of an optical beam propagating through turbulence. The filter fimction 
turns out to be the short-term modulation transfer function (MTF) . It is 
assumed that the log-amplitude covariance of a spherical wave for strong 
turbulence is given by the first-order log-amplitude covariance whose spatial 
frequency spectrum has been modified by multiplication with the spectral filter 
function (or the short-term MTF). The log-amplitude covariance thus calculated 
agrees well with the experiments and is shown in Figure 31. 
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From this covariance function, Clifford and Yura [60] obtained the 
following asymptotic behavior for the log-amplitude variance: 



0.36 (a^) 


-V, 


5 / 

s + 0.42 g: ^ 


( 50 ) 


where o' is a constant found from experiments. For a = 0.7, which is obtained 

from the strong scintillation data, tends to 0.23 for large values of 

X T 


The perturbation theory of Tatarski [32,34] implicitly assumes that the 
initial coherence is maintained along the propagation path and eddies having 
scale sizes of the order of the Fresnel length are most effective in producing 
scintillation. However, as the wave travels through the medium, the transverse 
spatial coherence decreases continuously and the wave becomes partially 
coherent. As was previously noted, the spreading of a beam of initial diameter 
d in a turbulent medium is characteristic of an aperture of diameter pq (the 
lateral coherence length) and not of diameter d. The radiator has become par- 
tially coherent in the presence of turbulence. Incorporating the loss of coherence 
of a radiator in the turbulent medium into the analysis, Yura [61] obtained the 
saturation of scintillations and the behavior of the amplitude covariance. He 
showed also that in the strong turbulence regime, the amplitude correlation 
length is equal to the lateral coherence length and increase. The turbulent 

eddies most effective in producing scintillations in the saturation regime are 
those that have scale sizes of the order of the lateral coherence length and not 
the Fresnel length. 


The asymptotic log-amplitude variance behavior of ( a^) ® has been 


predicted also by Fante [62] and Gochelashvily and Shishov [63] from approxi- 
mate solutions to the equation satisfied by the fourth moment of the field. 


The physical theories [59,60,61] and the approximate solutions [62,63] 
have been able to predict the observed behavior of the covariance function, 
namely, a rapid fall -off followed by a slow decay resulting in a long tail. These 
curves indicate that there are small-scale and large-scale structures in the 
scintillation pattern for strong turbulence. The correlation in strong turbulence 
is over a transverse separation of the order of 


Po 


( AL )'^2 


(51) 
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As increases, the lateral coherence length and the amplitude correlation 

length decrease. The diameter of the detector must be smaller than pg for 
point detector performance in strong turbulence. 


The aperture averaging which produced a reduction in scintillation noise 
became less effective for a given aperture of the receiver as the turbulence 
increased. This effect is due to the presence of large-scale scintillations at 
the receiver as increased. 


Another quantity of interest is the temporal frequency spectrum of inten- 
sity fluctuations observed at a fixed point. The width of the frequency spectrum 

1 / 

is of the order of V/(AL) ^ for weak turbulence where V is the velocity com- 
ponent normal to the beam. This is typically 10 to 100 Hz. For the case of 
strong turbulence, the width of the frequency spectrum is of the order V/p^ 


which is approximately 100 to 1000 Hz. In fact, the peak of the spectrum occurs 

1 / 

at approximately V/(\L) ^ for < 1 and at approximately V/p^ for > 1. 


In addition to the shifting of the peak of the spectrum to higher frequencies, 
there is a broadening of the spectrum and the maximum value of the spectrum 
decreases for increasing ct^. These features of the spectrum have been quali- 
tatively observed in the experiments. Finally, to obtain all of the temporal 
information of the scintillation in the strong turbulence regime, it may be 
necessary to use an electronic bandwidth of the order of V/p^. Otherwise, the 

use of a smaller electronic bandwidth may result in an apparent decrease of the 
log-amplitude variance for increasing values of a^. 


D. Effect of Turbulence on Heterodyne Systems 

In a heterodyne receiver, the incoming wave which is usually weak is 
combined with a reference wave on the photodetector surface. Coherent detec- 
tion results from the square-law response of the photodetector to the incident 
radiation. The squaring operation of the photodetector produces a current at 
the intermediate frequency which contains the signal -modulation information. 

The heterodyne detection of an atmospherically distorted wavefront has 
been the subject of several studies [64,65J. The various effects of atmospheric 
turbulence on the laser beam are related to the image behavior in a heterodyne 
system. Comparison of the expression for the optical power density due to a 


63 



laser transmitter at a point some distance L away from the transmitter and the 
expression for the S/N ratio of a heterodyne receiver detecting a signal from a 
point source at a distance L shows that they have identical functional dependence 
on the turbulence parameters. Thus, the telescoj e performance measured in 
terms of the coherence diameter remains the same whether it is used as part 
of a transmitter or as part of a heterodyne receiver for propagation in a turbu- 
lent atmosphere [66]. In simple terms, the reciprocity between the transmitter 
and receiver indicates that the average target illumination is related to the 
power received by the heterodyne receiver. Furthermore, beam wander and 
beam spread in the target illumination system are related respectively to image 
dancing and spread in the heterodyne receiver. Thus, the beam-image reciproc- 
ity implies the familiar effects of turbulence on the beam propagation in the 
analysis of the heterodyne detection. 

Fried [64] considered the amplitude and phase fluctuations and found 
that for an atmospherically distorted plane wave and a plane-wave local oscilla- 
tor there is an upper limit to the achievable S/N no matter how large the detector 
aperture. The aperture diameter for the near-peak value of S/N is of the order 
of the lateral coherence diameter. This aperture size has the virtue of avoiding 
large signal power fluctuations and variance which would result for larger collec- 
tion apertures. Moreland and Collins [65] tried to determine the possibility of 
increasing S/N by properly shaping the local oscillator beam and came to the 
conclusion that, though the optimum location oscillator shape differs significantly 
from that of a plane wave, the increase of average S/N over the plane-wave case 
is negligible. 


E. Signal-to-Noise Ratio Loss 

The wave function at a distance L from the transmitting telescope of 
radius R may be written by the Fresnel diffraction theory in the absence of 
turbulence along the path as [67,68,69] 


!/'(x,y,L) 



[exp i(KL - cot) ] 


+C0 p 

/ / exp i^ 

-CO L 


-ZL «ii2 
Af 



dx" dx” 


(51) 
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where 


= number of photons transmitted per second 
r” = vector in the plane of the transmitting lens 

r = vector in the plane at distance L from the transmitter and normal 
to the optic axis 

X = wavelength of laser radiation 
oj = angular frequency 
t = time. 


IS 


Equation (51) may be integrated by completing the square and the result 
^(x.y.L) = (^) + /(I . jl)) exp[^i(KL -ut+rt 

(gy)’’ 1 

r2[i + /g2(i - ^f] J 


where 


P = 


7TR^ 

XL 


C = T 


and 




( 52 ) 


( 53 ) 


( 54 ) 


( 55 ) 
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The effect of the atmospheric turbulence along the path of propagation is 
to produce random variations of phase and intensity. These random variations 
are introduced into the analysis by multiplying equation (51) with 


[X( r, r”) + i S( r, r”) ] 


(56) 


X( r , r”) is the perturbation of logarithm of amplitude, equation (41) , and 
S(T, 7”) is the perturbation of phase between the point 7" at the transmitter 
and the range point 7. Thus, the outgoing wave function at a distance L with 
the effect of turbulence is given by 


!/'(x,y,L) 



exp [i(KL 



( 7 - r ")2 


- i r"^ + i S( r, r”) + X( r, r") 


- dx” dy” 


(57) 


The transmitted wave is assumed to be backscattered by a diffuse target 
of scattering particles. Considering a single particle, the incoming wave is 
scattered by the particle as a spherical wave, and the scattered wave function 
at a point 7’ in the plane of the receiving lens for a coaxial system is given by 


i/'gC r’) 


r nV2 


(*) 


ip{ r) exp 


- (KL + Awt) - S( r, r') Ux(r, r') 


!-ir— 7’2 - 

( [Af AL 


(r - r’) 


T\ 2 


(58) 


where 


ct' = backscattering coefficient of the particle 

Aco = 2ojV/c = Doppler shift due to the velocity component V of the 
moving particle along the optic axis 

X( r , r') = perturbation of the log-amplitude 

S( r , r’) = perturbation of phase, between the range point r and the 
point 7’ at the receiver. 
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The local oscillator reference wave function which is added to the 
received wave function is assumed a plane wave with gaussian amplitude varia- 
tion which may be written as 




ref 


a A exp 



(59) 


where a is a constant and the arbitrary phase factor is assumed to be zero. 
The total signal photo-current is given by the real part of 


i 

s 


27) / If 


dr- 


(60) 


where W(r’ ) is the apodization or transmission function for the receiver aper- 
ture and is assumed equal to unity over the width of the reference beam. 

Substituting equations (57), (58), and (59) into equation (60) gives the 
signal current as the real part of 


i 

s 


\l 2(j^o r r 

ii 



+ i 


jAcot + 0( r*’ 


r’ 


r) 


where 


+ S( r, r") + S( r, r’) 


+ X( r, r”) + X( r, r’) 


dr' dr" 


(61) 


(/) ( r, r' , r”) = 2KL + 


XL 


( r - r' ) ^ + ( r - r" )^ 


^ (T' ^ + r"^ ) 


(62) 


The photo-current due to a single scatterer may be written as 


67 
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J/exp 
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'. r') 


+ X( r, r”) 


+ S(T, 7”) 


cos 


Awt + <p ("r”, r’ , r’) + S( V, r’ ) 


dr’ dr” 


( 63 ) 


The signal power is proportional to i*^ and 
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.9 2 ( 7 ’ r r r r 

^s = JjJJ 


j Jit 2 ^ t2 ^ 1 1 t 2 j»T TT t2 


+ X{r,^') + X(T, r”) + X( r, r’”) + X( r, r””) 




cos Acot 
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+ X( r, r') 


+ X( r, r”) + X( r, r' ”) + X( r, r’"') 


I cos j^2Acijt 


+ <p( r”, T’ , r) + ip( r'"', r' ”, r) + S( r, r' ) + S( r, r”) 

+ S(7, r’ ”) + S(T, r’"')J + cos r”, r' , r) + S( r, r' ) 


+ S(r, r”) -</)( r’"'+7'”' + 7) - S( r, r”') 
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If the signal power is averaged over a period longer than l/(2Aai) , then 


i^ = Re 
s 


r rr r r' 

^lW L" “ 


2 + j.tt2 ^ j,i ft2 ^ j,Mti2 


R^ 


+ X( r, r’ ) 


+ X( r, r”) + X( r, r' ”) + X( r, r'”') exp 


<^(r”, r', r) 


-0(7"”, r'", r) +S(r, r') + S(7, r") - S( r, r’") 


-S(7, r"") 


>dr' dr" dr’" dr"" 
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where Re denotes the real part. The point r' " lies on the receiver aperture and 
7 "" lies on the transmitting aperture. 

The expected value of the signal power is obtained by taking the ensemble 

average of i^ as given by 
s 



I in 


2 ^ j,tt2 J.T tt2 ^ j,itti2 

R^ 


+ i< 0( r", 7 , r) - 0( r""i 7 ", r) 


exp 


X(7, 7 ) 


+ X( r, r’") + X(7, r") + X( r, 7"") [ + i{s( r, r’)+ S(7, r") 


- S( r, r’") - S( r. 
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For this report it is assumed that the outgoing and incoming paths are 
statistically independent. (A more accurate treatment assuming gaussian dis- 
tributed random variables will be described in a later report. ) This assumption 
allows the following; 
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exp 


X(T, T’) + X(T, r"') } + i\S( r, r’) - S( r, r”') 


+ X( r, ■?’) + X(7, T””) + i S( r, r”) 



-S(r, P’*’) 


+ D( |r” - 


= exp 


D(|r’ - r'”|) 
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where D(r) is the wave structure function. Kolmogorov theory of isotropic 
turbulence gives a 5/3-power variation for the wave structure function: 


D(r) = 



( 68 ) 


where r is proportional to the coherence radius and has the value 

SL 


r = 0.088 A C , for plane wave 

a n ’ 


6 / - 3 / “ 6 / 

r = 0.1586 ^ , for spherical wave 

a n 


(69) 
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It will be appropriate to use r for a spherical wave for the wave originating at 
the range point. 

To simplify the calculations, the structure function is approximated as 


D(r) = 



( 71 ) 


and the plane-wave value of r is used for both ways to give a conservative 
estimate. 
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The expected value of becomes 
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Performing the integrations approximately, the expected value of the 
S/N ratio for a single scatterer is obtained as [67] 



nA^cr'TrR^T 

+ a^ + 


exp 



r+aM^^l - 0 


T 


(73) 


where a = (R/r ) and t is the pulse duration. Integrating equation (73) over a 

plane normal to the axis at the range L, the S/N ratio due to a plane target 
(dropping the angle brackets) is obtained; 


S_ _ Tjcr ’ tR^ A^ 

^ 4 L^(l + a^) + (f - L)2 


(74) 


In the absence of the turbulence effects, a = 0 and thus the S/N is reduced by 
the ratio 


P 


L^(l + a^) + (f - 
[f - L'' 

(t L; 


+ 


(75) 
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due to turbulence. If the system is focused to infinity, then 


P 


L^(l + s?) + 




l} + 


7t2r4 


( 76 ) 


Thus, the calculation of the loss of power due to turbulence reduces to 

the calculation of a = R/r and P. 

a 

The wave structure function for a plane wave is given by 


D(r) 



/ 

Path 


dL 

n 


(77) 


For a propagation path which is inclined at an angle 0' with respect to the 
horizontal 


h = H + L sin 0' 


(78) 


where h is any altitude above sea level and H is the altitude of the scattering 
surface above sea level. For flights at Edwards AFB, the dry lake is approxi- 
mately 600 m above sea level. The strength of turbulence is assumed to vary as 

< = (79) 


2 2 
where C _ is a constant, h' = h - H, and = 3200 m. The value of C 

nO 0 nO 

obtained from earlier measurements of Huffnagel and Stanley [70] , is 4.2 x 10“*'*. 
We use H = 600 m and assume that equation (79) is valid from h = 10 m. Sub- 
stituting equation (79) into equation (77), the structure function may be written 
as 
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( 80 ) 


D(r) = 3.38X 

' nO 


cosec 0' 




where l(i) is the incomplete gamma function given by 


r 



^ - 1 / 

/ u 'ie-“du , 
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h - H 
ho 


(81) 


and ij and are dimensionless lower and upper altitudes. From equation (80) 


r = 
a 




-3/5 


(82) 


For the flight tests, 0=7® and A = 10.6 pm. Using these values, r 


becomes 


r = 5. 8 X 10 
a 


-10 


'=no{«V -«V 


-3/5 


(83) 


The optics diameter is taken as 25.4 cm and this gives R = 0. 127 m; therefore. 


a = — = 2. 19 X 10^ r I(i ) - I(i )1 

r 1_ ^ 2^ ' l^'J nO 


(84) 


The incomplete gamma functions may be obtained from the Pearson Tables [71] . 

The calculations are made for the flight tests between 600 m and 4000 m 
2 -14 

altitude and for C „ = 4. 2 x 10 . The calculated S/N values shown in Figure 

nO 

32 also include the transmission loss. 
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Figure 32. Comparison of measured and theoretical S/N values. 


The loss due to the atmospheric turbulence is smaller than that due to 
absorption at 10.6 pm wavelength. As observed earlier, the absorption loss will 
be least at the DF laser wavelength, 3.9 pm, but the turbulence loss may be 
higher due to the smaller lateral coherence diameter. 
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A comparison of the losses due to absorption and turbulence over a 20 km 
horizontal path at altitudes of 2. 5 km and 5 km is given in Table 6 at CO 2 and DF 
laser wavelengths using the same optic size. Due to the smaller coherence 
diameter at DF wavelength, the optics diameter should be less than that at CO 2 
wavelength. By using the same optics size for both, the turbulence losses have 
been overestimated for the DF laser. Even so, the total loss of DF system is 
approximately 15 dB less than for CO 2 system at 2. 5 km. At 5 km altitude, the 
magnitude of the losses has dropped and there is an improvement of only 5 dB 
using the DF system. For altitudes below 5 km, the DF system is seen to give 
less total loss due to the atmosphere. 

TABLE 6. COMPARISON OF LOSSES AT CO 2 AND DF LASER 
WAVELENGTHS OVER A 20 km HORIZONTAL PATH 
AT 2. 5 km AND 5 km ALTITUDES 


Altitude 

Laser 

Absorption 

Loss 

(dB) 

Turbulence 

Loss 

(dB) 

Total Loss 
(dB) 

2. 5 km 

COo 

24 

2.6 

26.6 






DF 

1. 75 

9.05 

10.8 


CM 

0 

0 

9. 6 

0.37 

10.0 

5 km 


0.46 



DF 

4.9 

5.4 


Since the atmospheric conditions are constantly changing, it is desirable 
to estimate the maximum and minimum attenuation of the CO 2 laser radiation 
that may be expected due to the atmospheric propagation. The upper limit on 
the attenuation is estimated using the AFCRL Midlatitude Summer Hazy Atmos- 
phere and an exponential variation with altitude of the turbulence strength with 


2 -14 2/3 

C . = 4. 2 X 10 m 
nO 


The lower limit on the attenuation is estimated using 


AFCRL Midlatitude Winter Clear Atmosphere and an exponential variation with 

2 -15 -2/3 

altitude of the turbulence strength with C^^ = 4.2 x 10 


m 


The results 


for a two-way horizontal path of 20 km range at each altitude up to 10 km altitude 
are shown in Figure 33. The difference in the upper and lower limits of attenua- 
tion is considerable near sea level mainly because of the differences in the tem- 
perature, humidity, and aerosol concentration. The breakup of the attenuation 
into the transmission part and turbulence part at each altitude is given in 
Table 7. 
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ATMOSPHERIC ATTENUATION (dB) 



Figure 33. Atmospheric attenuation at 10.6 /xm. 



TABLE 7. SIGNAL LOSS IN ATMOSPHERIC PROPAGATION AT 
10.6fxm DUE TO ABSORPTION, SCATTERING, AND 
COHERENCE DEGRADATION FOR A TWO-WAY 
HORIZONTAL PATH OF 20 km AT 
EACH ALTITUDE 



AFCRL Midlatitude Summer Hazy Atmosphere 

2 -14 

and C „ = 4. 2 X 10 
nO 


Absorption and 

Turbulence 



Scattering Loss 

Loss 

Total Loss 

Altitude 

(dB) 

(dB) 

(dB) 

10,0 m 

68.0 

13.0 

81.0 

2. 5 km 

24.6 

2.6 

27.2 

5. 0 km 

9.6 

1.05 

10.65 

7. 5 km 

5. 8 

0.34 

6. 14 

10. 0 km 

3. 53 

0. 12 

3.65 


AFCRL Midlatitude Winter Clear Atmosphere 

2f -15 

and C ^ = 4. 2 X 10 
nO 

10 m 

17.3 

3.4 

20.7 

2 km 

10.5 

0.3 

10.8 

4 km 

6. 5 

— 

6. 5 

7 km 

3. 5 

— 

3. 5 

10 km 

1.8 

— 

1. 8 
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IV. AEROSOL EFFECTS 


A significant amount of laser energy in the 10 ;xm region is backscattered 
by the atmospheric aerosols. The exact amount is dependent upon aerosol num- 
ber density and size distribution. Although the size distribution of aerosols 
tends to remain unchanged with altitude, the number density decreases rapidly. 
This effect tends to counteract the advantage obtained with the decrease in 
attenuation with altitude. 

The effects of the aerosols appear both in the backscatter and extinction 
coefficients. The backscatter coefficient, /S, influences the radar cross-section 
for a volume distribution of scatterers, Uq, in the following manner [72]: 


where c is the speed of light and t is the duration of the transmitted pulse. The 
extinction coefficient, K, may in turn be separated into molecular and aerosol 
components as follows: 


K 




( 86 ) 


where K is the molecular extinction coefficient and y is the contribution due to 
m ' 

the aerosols. 

These effects may be included in the S/N ratio given by equation ( 74) . 

In the present expression, the effects of turbulence are ignored, the system is 
focused at infinity, and a’ is replaced by Uq, for the volume of scatters; addi- 
tionally, the aperture radius, R, from equation (74) is replaced by R = D/2 


S_ 

N 


7?Ect/37T 



exp (-2KL) 


(87) 
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A. Calculation of the Aerosol Backscatter and 
Extinction Coefficients 


The backscatter coefficient, )3, and the extinction coefficient, y, of 
aerosol particles when illuminated by radiation of wavelength X are given by the 
following integrals; 


CO 

/3 = J o-(a,X, m) n(a) da 
0 

CO 


The term a (a, X,m) is the backscatter cross-section of a particle of radius, a, 

and refractive index, m. The term Q„ ^(a, X, m) is the extinction cross-section 

Ext' 

of a particle of radius, a, and refractive index, m. The remaining term, n(a) da 
is the number of particles per unit volume between radii a and a + da. 

A Mie scatter program [73] has been used to calculate the scattering and 
extinction cross-section at each value of the radius, a, of the aerosol particles. 
The refractive index, m, is generally a complex number. In this case we use 
m = 1. 72 - i 0. 07 and m = 1. 6 - i 0. 12 in the calculations. 

There are a number of models for the particle size distribution n(a) . 

We computed /3 and y on UNIVAC 1108 for Haze L, Haze M, Log-normal, and 
Modified Haze C distributions and the results are presented in Table 8. All of 
these distributions are normalized such that the total number of particles per 
unit volume is 100. That is, we assumed in our calculations that 


J n(a) da = 100 . (88) 

0 


The details of the computer program are presented in Reference 74. 
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TABLE 8. BACKSCATTER AND EXTINCTION COEFFICIENTS AT 
10. 59 /urn FOR N = 100 particles/cm^ FOR 
SEVERAL DISTRIBUTIONS 


Distribution 

Refractive 

Index 

/3 

7 

Backscatter 

Haze L 

1.62 - 0.12 i 

0.1968 X 10“^ 

0. 1526 X 10“2 

Haze M 

1. 62 - 0. 12 i 

0.1905 X 10"3 

0. 1151 X 10“^ 

Log-Normal 

1.62 - 0.12 i 

0.2339 X 10“'' 

0.4533 X 10“® 

Haze C 

1.62 - 0.12 i 

0.3408X 10“® 

0.4233 X 10"® 

Extinction 

Haze L 

1. 76 - 0.07 i 

0.27 X 10“'* 

0.1053 X 10"® 

Haze M 

1. 76 - 0.07 i 

0.265 X 10"3 

0.1164 X 10"^ 

Log-Normal 

1.76 - 0.07 i 

0.7952 X 10"^ 

0.4342 X 10"® 

Haze C 

1.76 - 0.07 i 

0.73 X 10"® 

0.4214 X lO"® 


In view of the normalization given in equation ( 88) , we may write 

(3 = N(3 and y = Ny (89) 


where N is the total number of particles per unit volume, /3 is the backscatter 
coefficient for N = 1, and y is the extinction coefficient for N = 1. The relations 
presented in equation (89) express the assumption that the size distribution of 
aerosol particles remains the same and only the number density of aerosol 
particles changes from place to place. These results can now be used to 
calculate the behavior of the S/N ratio at ground level. 
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B. Performance of the MSFC Pulsed Laser System at Ground Level 

The dependence of the S/N ratio on aerosol particles is contained in the 
expression /3 exp (-2KL) which, using equation (89), maybe written as 

A = /3exp(-2KL) = N^exp[ -2(K^ + Ny) L ] . (90) 


If there are more aerosols in the atmosphere, N i . large and we get more back- 
scattered energy. However, aerosol extinction increases with N exponentially. 
Thus, it is not advantageous to have a very hazy atmosphere. The optimum 
value of the number density of aerosols for a given range may be found by 
differentiating equation (90) with respect to N and setting it equal to zero. The 
result is 


opt 


1 

29L 


(91) 


For modified Haze C distribution, we have y= 0.42 x 10~®/km. For example, 
we get = 1. 185 X 10^/cm® for a range of 10 km. Using the altitude distribu- 
tion of aerosol number densities given by the AFCRL model [75] (Taole 9) , the 
values of the range at which the S/N ratio reaches a maximum value due to the 
aerosols in the lower atmosphere are presented in Table 10. It may be observed 
that the optimum range at ground level in hazy atmosphere is limited to approxi- 
mately 9 km due to the large number of aerosols. As the altitude increases, 
aerosols do not pose any limitations on the S/N ratio due to the rapid decrease 
in their number density. 

We now consider the maximum possible ground range. Examination of 
equation (87) reveals that there is an optimum optics diameter for each range 
which maximizes the S/N ratio. Figure 34 illustrates the dependence of the 
S/N ratio on the diameter. Differentiation of equation (87) with respect to the 
diameter D and setting equal to zero results in the following relation: 


L 



4A 


(92) 
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TABLE 9. AFCRL MODEL OF AEROSOL NUMBER DENSITY 
VARIATION WITH ALTITUDE 


Altitude 

(km) 

Clear Atmosphere 
(particles/cm®) 

Hazy Atmosphere 
(particles/cm®) 

0 

2.828 X 10^ 

1.378 X 10^ 

1 

1.244 X 10^ 

5.03 X 10® 

2 

5.371 X 1Q2 

1.844 X 10^ 

3 

2.256 X 10^ 

6. 731 X 10^ 

4 

1.192 X 10^ 

2.453 X 10^ 

5 

8.987 X 10 

8.987 X 10 


TABLE 10. THE RANGE AT WHICH S/N REACHES A MAXIMUM 
DUE TO AEROSOLS HAVING HAZE C SIZE DISTRIBUTION 


Altitude 

(km) 

Clear Atmosphere 
(L - km) 

Hazy Atmosphere 
(L - km) 

0 

41.9 

8.6 

1 

95.3 

23.6 

2 

220.6 

64. 3 

3 

525 

186 

4 

994 

483 

5 

1319 

1319 
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Equation (90) is presented in numerical form in Figure 35. The corresponding 
S/N ratio is given by 



t]Ec 

4hu 



-2KL 

e 


(93) 


(S/N) is the maximum value of the S/N ratio obtained by using the optimum 
' ''max 

diameter of the optics. Any laser system with fixed optics cannot have a S/N 
ratio value higher thanthat of equation (93) . Thus, it forms an upper limit to 
the S/N ratio. 

Using equation (91), we obtain the ground range as follows, assuming 
that the minimum detectable S/N ratio is 10 dB. Taking the logarithm of both 
sides of the equation and with 10 dB, we obtain the relation between 

the extinction coefficient K in units of dB Am and the maximum range L : 

m 


K 


5 

L 

m 


log 


t]Ec A/3t 
4hi/L 

m 


dB/km 


(94) 


Setting the minimum detectable S/N equal to zero, we obtain: 



t]Ec X I3t 
4hfL 

m 


(95) 


The situation of equation (95) is not practical but is presented here for the 
purpose of comparison. 

At present, for the MSFC Laser System t^ = 0.05, E=10 mJ, and t = 

2 X 10“® s. The backscatter coefficient is left as a parameter varying between 
10“^ to 10“Vm. Based on these parameters, the maximum ground range 
obtained from equation (94) is presented in Figures 36 and 37. 

It may be observed from Figure 37 that the ground range is quite sensi- 
tive to the attenuation. A major component of the atmospheric attenuation is the 
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molecular absorption. The total atmospheric attenuation consists of the molecu- 
lar absorption, aerosol absorption, and scattering. Aerosol absorption and 
scattering amount to approximately 10 percent of the total attenuation at sea 
level and drops off rapidly with attitude. Therefore, molecular absorption alone 
gives a fairly good picture of the attenuation being experienced by the laser beam. 
Molecular absorption is extremely sensitive to the temperature and relative 
humidity. Using the results calculated in Section n, the molecular absorption 
is presented in Figure 38 as a function of temperature and relative humidity. 

This figure is useful in determining the attenuation rapidly when the temperature 
and humidity are known at the time of field test. 

The contribution of atmospheric turbulence discussed in Section III is 
usually in the range of 15 percent (or less) of the total atmospheric attenuation. 

In the calculation of the ground range, the turbulence effect may be included, 
but its effect on the ground range is not noticeable. 
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V. CONCLUDING REMARKS 


This report has discussed some of the various losses associated with the 
transmission of CO 2 laser radiation through the atmosphere, in particular 
molecular absorption, turbulence, and aerosol absorption and scattering. 
Additionally, calculations have been performed for selected data from test 
flights of the MSFC Pulsed Laser Doppler System. 

Specifically, the atmospheric transmission loss was calculated for five 
P lines (P16, P18, P20, P22, and P24) of the CO 2 laser using the January 1973 
fli^t test conditions. It was found that the transmission loss was approximately 
7 percent higher for P16 lines and 10 percent lower for P24 lines compared to 
the P20 line. Thus, the variation of the transmission loss is not very significant 
around the P20 line which is the usual range of efficient operation of the CO 2 
laser. The P40 line or RO line have better transmission, but it may not be 
possible to operate on these lines. Also, the attenuation of the CO 2 laser radia- 
tion due to absorption, scattering, and turbulence of the atmosphere is estimated 
for the January 1973 flight test conditions. The theoretical values of S/N, 
including the atmospheric attenuation, compare favorably with the flight test 
results. 

Molecular absorption was calculated for carbon dioxide, water vapor, 
and nitrous oxide for varying pressures and temperatures as a function of 
several P lines. 

The effects of atmospheric turbulence are analyzed and included in the 
expression for the S/N ratio along with aerosol effects relative to backscatter 
and extinction. 

The upper and lower limits on the atmospheric attenuation at each alti- 
tude up to 10 km are provided in graph form for a two-way horizontal path of 
20 km range. These calculations are based on the currently accepted models 
of atmospheric pressure, temperature, humidity, aerosol concentration, and 
turbulence strength. These results are expected to be useful in the interpreta- 
tion of the future test results at several altitudes using the Pulsed CO 2 Laser 
Doppler System. 
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